THE INTEGRALS IN GRADSHTEYN AND RYZHIK. 
PART 7: ELEMENTARY EXAMPLES 
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Abstract. The tabic of Gradshteyn and Ryzhik contains some elementary 
integrals. Some of them are derived here. 



1. Introduction 

Elementary mathematics leave the impression that there is marked difference 
between the two branches of calculus. Differentiation is a subject that is systematic: 
every evaluation is a consequence of a number of rules and some basic examples. 
However, integration is a mixture of art and science. The successful evaluation 
of an integral depends on the right approach, the right change of variables or a 
patient search in a table of integrals. In fact, the theory of indefinite integrals of 
elementary functions is complete OH]. Risch's algorithm determines whether a 
given function has an antiderivative within a given class of functions. However, 
the theory of definite integrals is far from complete and there is no general theory 
available. The level of complexity in the evaluation of a definite integral is hard to 
predict as can be seen in 



(1.1) 



/ da; = 1, / e"^' dx = and / e"^" dx (|) 



The first integrand has an elementary primitive, the second one is the classical 
Gaussian integral and the evaluation of the third requires Euler's gamma function 
defined by 

(1.2) r(a) = / x'^-'^e^'^ dx. 



The table of integrals [S] contains a large variety of integrals. This paper contin- 
ues the work initiated in [U [7l [H [9l [TOl [11] with the objective of providing proofs 
and context of all the formulas in the table [S] . Some of them are truly elementary. 
In this paper we present a derivation of a small number of them. 

2. A SIMPLE EXAMPLE 

The first evaluation considered here is that of 3.249.6: 
(2.1) / {l-y^)P-Ux ^ 







p{p + 1) 
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The evaluation is completely elementary. The change of variables y = 1 — ^Jx 
produces 

(2.2) I^~2 f\pdy + 2 f\p-Uy, 

Jo JQ 

and each of these integrals can be evaluated directly to produce the result. 
This example can be generalized to consider 



(2.3) 



I{a) = ^il-x^f-Ux. 
Jo 



The change of variables t = .t" produces 

(2.4) I{a) = fl-i r t^^^-\l - t)P-^ dt. 







This integral appears as 3.251.1 and it can be evaluated in terms of the beta 
function 

(2.5) B{a,b)^ [ x^'-^l-xf-^dx, 



(2.6) I{a) ^ a-^B {p,a-^) . 

The reader will find in [TT] details about this evaluation. 

A further generalization is provided in the next lemma. 
Lemma 2.1. Let n G N, a, 6, c G M with be > 0. Define u = ac~ b^. Then 

Proof. Substitute y — b + c-^/x and expand the new term (y — 6)". □ 

3. A GENERALIZATION OF AN ALGEBRAIC EXAMPLE 

The evaluation 



/oo 



dx TT 



(1 + x^) V4 + 3x2 3 
appears as 3.248.4 in [5]. We consider here the generalization 

r°° dx 

(3.2) q{a,b) 
We assume that a, b > 0. 



c^) \/b + ax"^ 



The change of variables x — Vbt/y/a yields 

dt 



/ — 

Jo a 



(3.3) q{a,b) „ , , 

where we have used the symmetry of the integrand to write it over [0,oo). The 
standard trigonometric change of variables t = tan if produces 

r^^^ cos ip dip 



(3.4) q{a,b)^2,/E 



a cos^ ip + b sin p 
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Finally, u — sin ip, produces 

(3.5) q{a,b) = 2V^f ^ . 

The evaluation of this integral is divided into three cases: 



Case 1. a = b. Then we simply get q{a, b) — Ij 



Case 2. a < b. The change of variables s — u\/b — aj y/a produces (6 — a)w^ = s^a, 
so that 

.O.N , ,N 2 r ds 2 _i 

(3.6) q{a,b) = ^=^ — = -=^= tan c, 

Vb - a Jo 1 + Vb - a 



with c = ^b — aj y/a. 

Case 3. a > b. Then we write 

(3.7) q{a,b) = 2y^ 



1 du 



a — {a — b)u 



2 ■ 



The change of variables u = ^/as/\/a — b yields 
(3.8) q{a,b)= ^ 



\/a-b Jo 1 - s 



where c = \/a — b/^/a. The partial fraction decomposition 
now produces 



(3.10) q(a, b) = In 



\/ a — b \ ^/a — \J a — b 

The special case in 3.248.4 corresponds to a = 3 and 6 = 4. The value of the 
integral is 2tan~^(l/\/3) = -I, as claimed. This generalization has been included 
as 3.248.6 in the latest edition of ^Sj. 

We now consider a generalization of this integral. The proof requires several 
elementary steps, given first for the convenience of the reader. 

Let a, 5 G M with a <b and rt e N. Introduce the notation 

(3.11) /-/„(a,fo) 

Then we have: 



(a + bf^Y Vl + t^ 



Lemma 3.1. The integral /„(a, 6) is given by 

(3.2) /..(M)^/' ",-:''-;''" . 

Jo (a + aw^)" 

with a = b — a. 

Proof. The change of variables v = t/^/T+W gives the result. □ 
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The identity 

(3.13) (1 - v'y-' = (1 - v'y + (1 - v'^-' i-{a + av 

produces 

a 



2\n-l _ n ..2^r^ , n „.2^r^-l J ^ , „.,2\ _ ^ 

a 



(3.14) I^{a,b)^^ I ^dv + - ^^dv. 

The evahiation of these integrals requieres an intermediate resuh, that is also of 
independent interest. 

Lemma 3.2. Assume z G M and n e N U {0}. Then 

dx _ 1 /2n\ / tan-i 2^2^'= 1 \ 

(3-5) (i + ^2^2)«+i y ~z 2fc(\^) (l + z2)fc j • 

Proof. Define 

^^•^^^ ^"^"^ i (i + ^i)"+i ^~zl (l + yV^- 

Take derivatives with respect to z on both sides of p.l6p . The outcome is a system 
of differential-difference equations 

az z z 



Z ' Z(l + z2)"+l' 

Solving for a purely recursive relation we obtain (after re-indexing n n — 1): 
(3.18) 2nFn{z) = {2n - l)F„_i(z) + ^ 



(l + z2)»' 

with the initial condition Fo{z) = ^ta.n^^z. This recursion is solved using the 
procedure described in Lemma 2.7 of Ij. This produces the stated expression for 

Fniz). ' □ 

The next required evaluation is that of the powers of a simple rational function. 

Lemma 3.3. Let a, b, c, d be real numbers such that cd > 0. Then 

ax^ + bY ^ a" 4a" [d rT,^ f^c- ad\'' fn\ f2k -2 

dx — \ W — tan \/c/dy 

k=l 



a" 


4a" 


c" 




4a" 


n , 



bc-ad\ ^ fn\ /2k - 2\ t4 2^^ d^ 



Aad J \k 

Proof. Start with the partial fraction expansion 
(3.19) G{x) a, be -ad 



4ad / \fc/\fc-l 

. k-l 



cx'^ + d c cd cx'^/d+1 

and expand G(a;)" by the binomial theorem. The result follows by using Lemma 
[321 □ 

The next result follows by combining the statements of the previous three lem- 
mas. 
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Theorem 3.4. Let a, 6 G M+ with a < b. Then 

dt 



/„+i(a, 6) 



(a + ^/TTW 



4. Some integrals involving the exponential function 
In [5] we find 3.310: 

(4.1) / e-P^dx = -, forp>0, 

Jo P 

that is probably the most elementary evaluation in the table. The example 3.311.1 

dx _ ln2 
Jo 1 + ~ V 

can also be evaluated in elementary terms. Observe first that the change of variables 
t = px, shows that (|4.2p is equivalent to the case p = 1: 

This can be evaluated by the change of variables u = that yields 

du 



(^•4) ' ' nil + uy 

and this can be integrated by partial fractions to produce the result. The param- 
eter in (|4.2|) is what we call fake, in the sense that the corresponding integral is 
independent of it. The advantage of such parameter is that it provides flexibility 
to a formula: differentiating (j4.2p with respect to p produces 

r°° xeP'^dx _ ln2 
/o (l + eP-)2 " 



(4.5) 

i-^ x^eP'^ieP"' - l)dx 2 In 2 



(4.6) 



(4.7) 



(1 + eP^)3 p3 ' 

°° x^eP'^ie^P'' - 4eP^ + 1)^2; _ 6 In 2 



(l + eP^)4 p4 

The general integral formula is obtained by differentiating (j4.2p n-times with respect 
to p to produce 

The pattern of the integrand is clear: 

/ 9 \ " 1 (— l)"a;"eP^ 

(4-9) (^g^j iq:^ = (n_eP-)"+i^"(^''")' 

where Pn is a polynomial of degree n — 1 . It follows that 

x"eP^P„(eP^)(ix n!ln2 



(4.10) 



(1 + eP^)"+i 



,11+1 
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The change of variables t = px shows that p is a fake parameter. The integral is 
equivalent to 

(4.11) r-:-^^"(flf ^n!ln2. 

The first few polynomials in the sequence are given by 

(4.12) P,{u) = 1, 
P2{u) = u-l, 
P3(m) = u2-4u+l, 

P4{u) = - llu^ + llu _ 1. 

Proposition 4.1. The polynomials Pn{u) satisfy the recurrence 

(4.13) Pn+l{u) = [nu - l)Pn{u) - U{1 + u)^Pn{u). 

au 

Proof. The result follows by expanding the relation 

^^^^^ (-l)"+ia;"+ieP^P„+i(eP"^) d (-l)"a;"eP"^ ^^(eP^) 



(l + eP^)"+2 dp \ (l + eP^)"+i 



□ 



Examining the first few values, we observe that 

(4.15) g„(u) := (-l)"P„(-w) 

is a polynomial with positive coefficients. This follows directly from the recurrence 

(4.16) Qn+iiu) = [nu + l)Qniu) + u{l ~ u)^Qn{u). 

au 

This comes directly from (|4.13p . The first few values are 

(4.17) Qi{u) ^ 1, 

(32 (m) = U+1, 

Qa{u) = + llu^ + 11m + 1. 

Writing 



(4.18) Qniu)^J2^i 

3=0 

the reader will easily verify the recurrence 

(4.19) Eo,n+l = Eo,n 

Ej,n+i = {n-j + l)£^j-i.„ + (j + l)Ej^n 

The numbers Ej^n are called Eulerian numbers. They appear in many situations. 
For example, they provide the coefficients in the series 

rn — 1 



(4-20) E^^'-' = 7T^^E%"-" 

k=l ^ ' n=0 

and satisfy the simpler recurrence 

(4.21) Ej^n = nE.j^i^n + j£^j\n-l, 
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that can be derived from (I4.19|) . These numbers have a combinatorial interpreta- 
tion: they count the number of permutations of {1, 2, . . . , n} having j permutation 
ascents. The exphcit formula 



'n + 1 



(4.22) ^.■." = E(-1)'( k jO' + l-^)"' 

can be checked from the recurrences. The reader will find more information about 
these numbers in [6] . 

5. A SIMPLE CHANGE OF VARIABLES 

The table [5] contains the example 3.195: 
^^^ ^ Jo (X + a)P+i - p(a - 1) ■ 

One must include the restrictions a > 0, a I, p =/= 0. The evaluation is elementary: 
let 

(5.2) u-ii^, 
^ ^ x + a 

to obtain 

1 



(5.3) / uP~Uu, 

a - 1 Jl/a 

that gives the stated value. The formula has been supplemented with the value 1 
for a = 1 and In a/ (a — 1) when p = in the last edition of 

Differentiating (|5.ip n times with respect to the parameter p produces 

k=l 



/o {x + a)P+'^"' \x + aj (a-l)p"+i 
Naturally, the integral above is just 



(5.4) / uP-^ ln"udw 

« - 1 Jl/a 

and its value can also be obtained by differentiation of (|5.3|) . 

The next result presents a generalization of (|5.ip : 
Lemma 5.1. Let a, b be free parameters and n G N. Then 

Jo (x + «)"+" I ^0 



k J b + k 



where B{n, b) is Euler's beta function. 



Proof. Use the change of variables u = {I + x)/{a + x), expand in series and then 
integrate term by term. □ 
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6. Another example 
The integral in 3.268.1 states that 

To compute it, and to avoid the singularity at x = 1, we write 

(6.2) / = lim ' ( — dx. 

^ ' ^-oio \l-x 1-xPj 

This evaluates as 

(6.3) / = lim - Ine + ln(l - (1 - e)P) = lim In ( ] = inp. 

7. Examples of recurrences 

Several definite integrals in [5] can be evaluated by producing a recurrence for 
them. For example, in 3.622.3 we find 

(7.1) I w"... = ,-ir^j-g<^j. 

To check this identity, define 

^7r/4 

(7.2) /„ = / tan^^xdx 



and integrate by parts to produce 
(7.3) /„ - ^ 



ri-l 



2n- 1 

From here we generate the first few values 

TT^ 7r_^7r_l ,^ 7r_ll 
/o = ^, A = -4+1,/. = 4 -1+3, and/3 = -- + l-- + -, 

and from here one can guess the formula (|7.ip . A proof by induction is easy using 



A similar argument produces 3.622.4: 
(7.4) / tan2"+ia;dx = ^^^ 1^2-^^ ^ 



k 

k=l 



To establish this, define 

(7.5) J„ = / t£in'^'^+^ X dx 

Jq 

and integrate by parts to produce 

(7.6) J„ = - J„_i + 

2n 

The value 

(7.7) Jo— t'Anxdx — —— 

Ja 2 

and the recurrence (|7.6p yield the formula. 
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8. A TRULY ELEMENTARY EXAMPLE 

The evaluation of 3.471.1 

is truly elementary: the change of variables t = —b/x gives the result. 

9. Combination of polynomials and exponentials 
Integration by parts produces 

(9.1) / a;"e''^da; = -a;"e''^ - - / x"-^e'''=. 
J a a J 

This appears as 2.321.1 in Introduce the notation 

(9.2) /„(a) := J x^'e"'-' dx 
so that (19.11) states that 



1 r? 

(9.3) /„(a) = -x"e"^--/„_i(a). 

This recurrence is now used to prove 

{-Ifx""-^ 



(9.4) /„(a) = 7i!e"^^ 



by an easy induction argument. This appears as 2.321.2. The case \ < n < A 
appear as 2.322.1,2.322.2,2.322.3,2.322.4, respectively. 

Integrating (|9.4p between and u produces 3.351.1: 
(9.5) / ^™e-''-dx= ^-e-'^^y 

V I J n+l J.\ ^n-fe+l 

U fc=0 

This sum can be written in terms of the incomplete gamma function. Details will 
appear in a future publication. Integrating (j9.4p from u to oo produces 

This appears as 3.351.2. 

The special case n = 1 of 3.351.1 appears as 3.351.7. The cases n = 2 and n = 3 
appear as 3.351.8 and 3.351.9, respectively. 

10. A PERFECT DERIVATIVE 

In section 4.212 we find a list of examples that can be evaluated in terms of the 
exponential integral function, defined by 

I"" p* di 

(10.1) Ei(x) • 



for x < and by the Cauchy principal value of (jlO.ip for a; > 0. An exception is 
4.212.7: 

In X da; e 
/i (1 +lna;)2 ~ 2 



I 111 Ju U,Uy L 

(10-2) / 7^^rTT2 = o-l- 
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This is an elementary integral: the change of variables t — 1 + In x yields 

it -I) 



(10.3) 

and to evaluate it observe that 
(10.4) 



/.I 

e 



t2 



'-e^dt 



f2 



d e* 
diT' 



The change of variables t = In a; in (|10.2p yields 
(10.5) 

This is 3.353.4 in [S]. 



^ te^dt _ e 
(l + i)2 " 2 ~ ■ 



The previous evaluation can be generalized by introducing a parameter. 

Lemma 10.1. Let a e M. Then 

In a; da; e 



(10.6) 



_ e 1 
H (a + lna;)"+i ^ (a + 1)" a"' 

Proof. Substitute t — a + liix and use 

c? e* t — a t 

10.7) = — -re ■ 

The case a = 1 corresponds to (|10.2p . 

11. Integrals involving quadratic polynomials 



□ 



There are several evaluation in [5' that involve quadratic polynomials. We as- 
sume, for reasons of convergence, that the discriminant d = b'^ — ac is strictly 
negative. 



We start with 



dx 



1 



: cot 



^ ^ lo ax'^ + 2bx + c y/ac - 6^ \y/ac - b"^ , 

This is evaluated by completing the square and a simple trigonometric substitution: 

dx 1 [°° du 



ax"^ + 2bx + c 



a V? — dja 



2 

dv 



Differentiating (jll.ip with respect to c produces 3.252.1: 



dx 



(-1) 



n—1 Q'n—1 



(ax2 + 2bx + c)" (n - 1)! 9c"-i 



cot 

HVVac- &2) 
vac— 1? 



We now produce a closcd-from expression for this integral. 



Lemma 11.1. Let n e N and u := 4(oc — b^)/ac. Assume cu > 0. Then we have 

the explicit evaluation 

(11.2) 

dx 2b /2n — 2\ \ \/acu , / 2fe \ ^^-4 



/q (aa;2 + 26a; + c)" a(cu)" \ rt — 1 



■ cot 
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Proof. The case n = 1 was described above: 
(11.3) /i(a,6,c):=/ „ , , = ^ ^ cof^ 



Q ax^ + 2bx + c ^ac -Ip- \ Vac - P , 

Now observe that h{a^,abc,b^) = h{l,b,l)/ac. To establish pi.2p . change the 
parameters sequentially as a i-^ a^; ci—>c^; & i— > abc. In the new format, both sides 
satisfy the differential-difference equation 

(11.4) -2nc{l~b')U+, ■ ^ 



oc 



The identity (|11.2[) is obtained by reversing the transformations of paramaters 
indicated above. □ 



Corollary 11.2. Using the notations of Lemma [11.11 we have 

> 3 ■ 

The integral 
(11.6) 



(11-5) E-pTV^V^"^" 

7 = 1 JKj) 



dx (2n- 3)!!7ra"-i 



_^ (0x2 + 26a; + c)" (2n - 2)!! (ac - 62)»-i/2 



reduces via m = a(a; + b/a)/y/ac— 6^ to Wallis' integral 

(2n-3)!! tt 



(11.7) / 

Jo 



/o (w2 + l)" (2n-2)!! 2' 

that appears as 3.249.1. The reader will find in [2] proofs of Wallis' integral. 
Observe that the evaluation of 3.252.2 is much simpler than the corresponding 
half-line example presented in Lemma 1 11. II 

The last example of this type is 3.252.3: 

dx (-2)" a" / 1 



/o (ax2 -f 26a; c)"+3/2 {2n + 1)U dc'' + b) 

A simple trigonometric substitution gives the case n = 0: 

dx \/a f°° du 



(aa:;2 -|- 26a; 4- c)3/2 ac - 62 Jb/^^ [u^ + 1)^/2 

1 



Vc(Vac+ 6) 

The general case follows by differentiating with respect to c and observing that 

We now provide a closed- form expression for (jll.Sp . 

Theorem 11.3. Let a, 6, c € R and rt e N. Define u = [ac — b'^)/Aac and assume 
cu > 0. Then 



dx _ (cu)-'' ( 1 ^ ^'^^^ 



(aa;2 26a; c)"+3/2 ^{^^){2n + I) \^/^ + b ac-b^^^\j 



U-' 
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Proof. Change parameters sequentially as a i—> a^; c i—> c^; b abc. Then, in the 
new format both sides satisfy the differential-difference equation 



2^„^ J. _ dfN b 



(11.8) - (27V(1 - 5^)c)/jv+ 

where N ~ n + ^. □ 

12. An elementary combination of exponentials and rational 

functions 

The table fS] contains two integrals belonging to the family 

dx 



(12.1) Tj / e^P^(e~^ - 1)" 

Indeed 3.411.19 gives Ti: 



/ e-P^(e~^-l)"— = -y(-l)* 



(12.2) / e-^'-(e~--l)"— = -^^(-l)^(^^Jln(p + n-fc), 
and 3.411.20 gives Ta: 

(12.3) r e~P-{e-^ - 1)"^ = ^(-1)^= (l) {p + n - k) ln(p + n - k), 

Jo X fe^o v^y 



The next result presents an explicit evaluation of T, 



Proposition 12.1. Let p be a free parameter, and let n, j S N with n -\- p > 0. 
Then 



(12.4) / e-P-(e"--l)"— = -i-j-^(-l)'=(p + n-A;)J"-Mn(p + n-A;). 

Proof. Start with the observation that 

(12.5) Tj = - j T,^i{p)dp + C. 

Therefore we need to describe the iterative integrals fj{p) — J fj-i{p) dp, with 
fa{p) = \n{p + a). This can be found in page 82 of [2] as 

(12.6) f,{p) = ^{p + ay Hp + a)-^{p + a)^ + C, 

with a = p + n — k and Hj = 1 + ^ + ■■■ + ^ is the harmonic number. To build 
back the functions Tj employ the fact that, for any polynomial Q{n, fc), 

(12.7) j2(-l)\-lf(l)Qin,k)^0. 



Consequently, 



(12.8) T,=C+^ — ^ y{-iyip + n-ky\n{p + n~k). 



■' k=0 



The last step is to check that C = 0. This follows directly from Tj as p ^ oo. 
The assertion is now validated. □ 
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13. An elementary logarithmic integral 
Entry 4.222.1 states that 



(13.1) / \n[^^^]dx = {a-h)i:. 







In order to establish this, we consider the finite integral 



(13.2) I{m) := I In ( ] dx 







and then let to — s- cx). 

Integration by parts gives 



ln(a^ + x^) dx = TOln(TO^ + a^)-2 



X dx 
0? + x"^ 

dx 



= m\n(rn + a ) — 2m + 2a 
= TOln(TO^ + a^) — 2to + 2atan^^ ^ — ^. 

Therefore 

Urn) = TOln (^-^ j + 2atan- (-) - 26tan- (-) . 

The limit of the logarithmic part is zero and the arctangent part gives (a — 6)7r as 
required. 

The generalization 

(13.3) /°°lnf^l±^) dx = (a-6)— ^ 

^ ^ Jo yb'+x'^J ^ ^sin(7r/s) 

can be established by elementary methods provided we assume the value 

/■°° dx TT 

(13.4) 



1 + x* s sin(7r/s) 

as given. This integral is evaluated in terms of Euler's beta function in [llj . Indeed, 
integration by parts gives 

dx 

(13.5) / \n{a' +x')dx^y\n{a' +if)-sy + sa' 

Jo Jo + 2/* 

and similarly for the ^-parameter. Combining these evaluations gives 

' fa' + x'\ ^ , fa'+y'\ , /"^ dx T dx 



b'^ + x'' J \ b'' + J Jo '^^ + Jo + 

Upon letting y — + oo, we observe that the logarithmic term vanishes and a scaling 
reduces the remaining integrals to l|13.4p . 
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NSF-DMS 0409968. 



14 



TEWODROS AMDEBERHAN AND VICTOR H. MOLL 



References 

[1] T. Amdeberhan, L. Medina, and V. Moll. The integrals in Gradshteyn and Ryzhik. Part 5: 

Some trigonometric integrals. Scientia, to appear. 
[2] G. Boros and V. Moll. Irresistible Integrals. Cambridge University Press, New York, 1st 

edition, 2004. 

[3] M. Bronstcin. Integration of elementary functions. University of California, Berkeley, Cali- 
fornia, 1987. 

[d] M. Bronstein. Symbolic Integration I. Transcendental functions, volume 1 of Algorithms and 

Computation in Mathematics. Springer- Verlag, 1997. 
[5] I. S. Gradshteyn and I. M. Ryzhik. Table of Integrals, Series, and Products. Edited by A. 

Jeffrey and D. Zwillinger. Academic Press, New York, 7th edition, 2007. 
[6] R. Graham, D. Knuth, and O. Patashnik. Concrete Mathematics. Addison Wesley, Boston, 

2nd edition, 1994. 

[7] V. Moll. The integrals in Gradshteyn and Ryzhik. Part 1: A family of logarithmic integrals. 

Scientia, 14:1-6, 2007. 

[8] V. Moll. The integrals in Gradshteyn and Ryzhik. Part 2: Elementary logarithmic integrals. 
Scientia, 14:7-15, 2007. 

[9] V. Moll. The integrals in Gradshteyn and Ryzhik. Part 3: Combinations of logarithms and 
exponentials. Scientia, 15, 2007. 
[10] V. Moll. The integrals in Gradshteyn and Ryzhik. Part 4: The gamma function. Scientia, 15, 
2007. 

[11] V. Moll. The integrals in Gradshteyn and Ryzhik. Part 6: The beta function. Scientia, to 
appear. 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: tamdeberhaiiamath.tulane.edu 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: vhrnSmath . tulane . edu 



